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Abstract 

We calculate electromagnetic pion form factors with an analytic model for 
as(Q^) which is infrared (IR) finite without invoking a "freezing" hypothesis. 
We show that for the asymptotic pion distribution amplitude, F^o^*^ agrees 
well with the data, whereas the IR-enhanced hard contribution to F^^ and the 
soft (nonfactorizing) part can jointly account for the data. 
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The issue of computing exclusive processes, like electromagnetic form factors, within 
QCD is of fundamental interest because it reveals the basic structure of hadrons. But in 
contrast to inclusive processes, there is much uncertainty about the applicability of pertur- 
bative QCD at laboratory momenta. It seems that the more interpretations this subject 
inspires and grounds the more unsettled it is. 

This paper attempts to bring together some crucial results about pion form factors and 
combine them with novel theoretical developments concerning the infrared (IR) regime of 
QCD, in order to benchmark the status of our current understanding. Such an approach 
appears attractive since the incorporation of nonperturbative power corrections in the per- 
turbative domain may improve both the IR insensitivity of exclusive observables and the 
self-consistency of calculations entrusted to perturbative QCD |l|J^,3,[4|,^,|6|.|7|,|8|j9|Jl0(| . 

To this end, the conventional representation of the running strong coupling "constant" 
is given up in favor of an analytic model, recently proposed by Shirkov and Solovtsov [|l|, 
which incorporates a single power correction to remove the Landau singularity. Bearing in 
mind that the definition of beyond two loops cannot be uniquely fixed, one may regard 
the ambiguity in the IR modification of the running coupling as resembling the freedom of 
adopting a particular non-IR- finite renormalization scheme ||1 1|| . 

The Shirkov-Solvtsov model employs Lehmann analyticity to bridge the regions of small 
and large momenta by changing the one-loop effective coupling to read 



Po 



+ 



ln(QVA2) A2-g 



(1) 



where A = Aqcd is the QCD scale parameter. This approach provides a nonperturbative 
regularization at low momenta and leads to a universal value of the coupling constant at 
zero momentum a^^\Q'^ = 0) = 47r//9o — 1.396 (for three flavors), defined only by group 
constants, i.e., avoiding the introduction of external parameters, like an (effective) gluon 
mass to "freeze" the coupling at low momentum scales. 

Note that this limiting value (i) does not depend on the scale parameter A - this being 
a consequence of renormalization group invariance - and (ii) extends to the two-loop order, 
and beyond, i.e., a^^HQ^ = 0) = ^^(Q^ = 0) = a,(Q2 = o). (In the following the 
bar is dropped.) Hence, in contrast to standard perturbation theory, the IR limit of the 
coupling constant is stable, i.e., does not depend on higher-order corrections, and is therefore 
universal. As a result, the running coupling also shows IR stability. This is tightly connected 
to the nonperturbative contribution oc exp(— 47r//?o) which ensures analytic behavior in the 
IR domain by eliminating the ghost pole at = A.qqy). 

At very low momentum values, say, below 1 GeV, Aqcd in this model deviates from 
that used in minimal subtraction schemes. However, since we are primarily interested in a 
region of momenta which is much larger than this scale, the role of this renormalization- 
scheme dependence is only marginal. In our investigation we use AQ(lg = 242 MeV which 

corresponds to Aq^^'' = 200 MeV. 

This analytic model for the strong running coupling is very suitable for calculations of 
exclusive amplitudes, mainly for two reasons: Firstly, it ensures IR safety of the factorized 
short- distance part without invoking the additional assumption of saturation of color forces 
by using a gluon mass - extensively used up to now in form-factor calculations (see, for 



example, [0). Furthermore, the Sudakov form factor |jT3[ does not have to serve as an 
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IR protector against as singularities. Hence, the extra constraint of using the maximum 
between the longitudinal and the transverse scale, as argument of Og, proposed in [14 



becomes superfluous. This is a serious advantage relative to previous analyses because now 
one is able to choose the unphysical constants [|15|, which parametrize different factorization 
and renormalization schemes, in such a way as to optimize calculated observables (for a 
more detailed discussion of this point, we refer to Second, and more important, the 

non- logarithmic term in Eq. (|I|) enters all anomalous dimensions, viz. the cusp anomalous 
dimension, which gives rise to the Sudakov form factor |[T6|JT7|JT8|JT9|j20[] , as well as the quark 
anomalous dimension which governs evolution. As a result, the suppression due to transverse 
momenta, intrinsic ||21| and those generated by radiative corrections |T^, is counteracted, 
and hence there is no reduction of the form-factor magnitude. 

We are going to show in this work that the enhancement effect originating from the power 
correction to the running coupling is enough for the asymptotic pion distribution amplitude 
to contribute (at leading order) to the spacelike electromagnetic form factor of the pion, 
^iriQ'^), cL hard part that can account for almost half of the form-factor magnitude relative 
to the existing data [0,^ . 

On the other hand, the transition form factor, F^o^.^(Q^), is only slightly changed, as 
compared to the result given in |2^, and matches the recent high-precision CLEO data ||2^ 
as good as the dipole fit. Also the older CELLO data at lower Q"^ are well reproduced 
(see below). 

In both cases, no adjustment of the theoretical predictions to the experimental data is 
involved. 

Therefore, there appears to be no need to reanimate endpoint-concentrated pion dis- 
tribution amplitudes, proposed by Chernyak and Zhitnitsky (CZ) |2^, in order to make 
contact with the experimental data - as recently attempted in [^. We are less enthusiastic 
about using such distribution amplitudes because of the following serious disadvantages: 
(i) It has been recently shown P^|2U| that distribution amplitudes of the CZ-type lead to a 
TTj transition form factor which significantly overestimates the CLEO data just mentioned. 
Our reasons for skepticism parallel the arguments given in |24] and will not be repeated here. 
The excellent agreement between theory (QCD) and measurement for this process, already 
at leading order, when the asymptotic pion wave function is used, cannot be overemhasized. 
Note in this context that the calculation of Cao et al. ||30| 
function smaller values of F^o^*^((5^ 



which predicts for the CZ wave 
than the data, uses for modeling the k±_ distribution 
in the pion an ansatz that strongly suppresses the endpoint region. Hence, in effect, their 
wave function, though claimed to be the CZ one, excludes this region and yields therefore a 
result even smaller than the one predicted by the asymptotic distribution amplitude in the 
range of where there are data. Independently, investigations based on QCD sum 

rules come to a comparably good description of F^o^*^(Q^) at not too low on the basis 
of local duality without presuming the asymptotic form of the pion distribution amplitude, 
but favoring again a shape close to that. 

(ii) Distribution amplitudes of the CZ-type yield a direct-overlap, i.e., soft contribution, 

to F^(g2 



33.34.35.24 



as 



which turns out to be of the same large order of magnitude 
that resulting from the convolution with the hard-scattering amplitude [^,^,^], or is even 
larger, at currently probed values. Inclusion of this contribution into the pion form factor 
leads eventually to a total result which overestimates the existing data considerably - even 
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allowing for some double counting of hard and soft transverse momenta near the transition 
region. Though the present quality of the high-Q^ data |22,23] on the spacelike pion form 
factor is quite poor, the trend seems to be indicative. 

(iii) The underlying QCD sum rules analysis of ||2^ suffers with respect to stability - 
as outlined by Radyushkin ||3^. As a result, the duality interval increases with moment 
order N, meaning that the (nonperturbative) condensate contributions grow with order 
relative to the perturbative term. A direct consequence of this is that the moments of 
the pion distribution amplitude for = 2, 4, extracted from these sum rules, are artificially 
enhanced. Such large moment values can only be reproduced by a double-humped endpoint- 
concentrated distribution amplitude and correspond to the basic assumption that vacuum 
field fluctuations have infinite size, or equivalently that vacuum quarks have exactly zero 
virtuality ^4 



(iv) The characteristic humps in the endpoint regions (x = 0, 1) are not generic, but 
merely the result of truncating the eigenfunctions expansion of Gegenbauer polynomials 
at polynomial order two while keeping the normalization fixed to unity. Including higher 
and higher order polynomials, the humps become less and less prominent and the central 
region (x = 1/2) gets enhanced. To this point, we mention that an independent QCD sum 
rules analysis gives the constraint (pj^ix = 1/2) = 1.2 ± 0.3, which is close to the value 
(pas{x = 1/2) = 3/2, and definitely violated by the CZ amplitude. 

Physically, the source of the endpoint enhancement of CZ-type distribution amplitudes 
can be understood as follows. If the vacuum quark virtuality is zero, an infinite number 
of such quanta can migrate from the vacuum to the pion state at zero energy cost. This 
happens exactly in the kinematic region x = Oora; = l — x = and leads to a strong 
enhancement of that region at the expense of depleting the amplitude for finding configura- 
tions in which the quark and the antiquark, or more precisely, the struck and the spectator 
partons, share almost equal momentum fractions around a; = 1/2. In the pion, configura- 
tions close to the kinematic boundary contain one leading parton, which picks up almost all 
of the injected momentum, and an infinite number of wee partons ~ l/x with no definite 
transverse positions relative to the electromagnetic probe, which constitute a soft cloud. In 
this regime, gluons have very small virtualities and therefore it is inconsistent to assume 
hard-gluon rescattering, i.e., the factorization of a short-distance part in the exclusive am- 
plitude becomes invalid. This region of momenta has to be treated separately on the basis 
of the Feynman mechanism just described, but a theoretical approach from first principles, 
though of paramount importance, is still lacking. On the other hand, if the vacuum virtual- 



ity is sizeable, say, of the order of Aqcd or even larger |3^, then an energy gap might exist 
that prohibits the diffusion of vacuum quarks into the pion state, and hence Feynman-type 
configurations are insulated from those for which hard-gluon exchange applies. This gap 
may be the result of nonlocal condensates p4| , ^^ , which have a finite fluctuation size, 
or alternatively being induced in the form of an effective quark mass acquired through the 
interaction with an instanton background |^2|,^. But the general result is the same: the 



shape of the pion distribution amplitude gets strongly enhanced in the central region and 
resembles closely the asymptotic one. 

In view of these drawbacks, a potentially good agreement between theoretical estimates 



employing CZ-type distribution amplitudes for the pion - as recently reported in |^ - and 
experimental measurements is entirely circumstantial. 
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In the present work, we are going to show that taking together the soft form-factor 



contribution due to the overlap of the initial and final pion wave functions and the hard, 
i.e., factorizing part, we can obtain a result in qualitative agreement with the existing data 
and complying with the power counting rules. Actually, including also NLO contributions 
to the hard-scattering amplitude (see, and previous references therein), Ft,{Q'^) gets 



additionally enhanced to account for approximately half of the form-factor magnitude |]1T 
modulo the large uncertainties of the existing data. 

However, we cannot and do not exclude that the true pion distribution amplitude may 
deviate from the asymptotic one, but this deviation should be within the margins allowed 
by the experimental data for the 7r7 transition form factor. Hence the true pion distribution 



amplitude may well be a hybrid of the type $t 



90%$as + 9%$cz + 1%CT. This 
mixing ensures a broader shape of the pion distribution amplitude, with the fourth-order, 
"Mexican hat"-like, Gegenbauer polynomial, being included in order to cancel the dip at 
X = 1/2. The shapes derived from instanton-based approaches are of this type. For 

such distribution amplitudes, evolution already at LO must be taken into account that tends 
to reduce the importance of the endpoint region leading to a decrease of the magnitude of 
the form factors towards the data p4| , ^ . On the other hand, for the asymptotic solution, 
evolution enters only at NLO and is a tiny effect which is ignored in the present exploratory 
investigation. 

The starting point of our analysis is the expression for the pion form factor in the 
transverse configuration space after employing factorization to separate a short-distance. 



i.e., hard-scattering part (where the the terminology of |]T5|,|T3],|T4| is adopted): 



FJQ 



dxdy 



J-oo (47r)2 



VT (y, P\ /iF, /iR) Th (x, b, Q/fin, QZ/ip) 



(2) 



Here P+ = Q/V2 = p-', = -{P'-Pf, and /ir = C2iQ and /ip = Ci/h are, respectively, 
the renormalization and factorization scales, with Ci, C2-\/2 = C^*^ (^ = x,x,y,y), the 
constants Ci, C2 being integration constants of order unity, so that (uncalculated) higher- 



order corrections are small ||T5|,|T8|. Finally, b is the variable conjugate to the transverse 
gluon momentum, and denotes the transverse distance between quark and antiquark. 

The hard-scattering amplitude Th is the amplitude for a quark and an antiquark to 
scatter coUinearly via a hard-gluon exchange with wavelengths limited by b, and is given in 
leading order by 



Th (x, y, b, Q/fiK, bQ) = 8CFar(/^R)^o iV^bQ) 



(3) 



In Eq. (|2|), describes the valence qq amplitude which includes gluonic radiative correc- 



tions uM as well as the primordial transverse size of the bound state [21 1: 



(x,6, Q,/iF,yUR) = exp 



f>^^ d/i 

-s (x, b, Q, hf, /iR) - s (x, b, Q, /iF, /Ur) - 2 / — 7q {g{fi)) 
X P(x,6,/iF) . (4) 
The pion distribution amplitude at the factorization point is approximately given by 
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V (x, 6, /ip = Ci/h) ^ (pT, (x, /ip = Ci/h) S(a;, 6) , 



(5) 



where 6) parametrizes the intrinsic transverse size of the pion (see below). In the 
coUinear approximation, one has 



(6) 



where /^r = 130.7 MeV and = 3. Integrating on both sides of this equation over x 
normalizes 0,^ to unity, i.e., Jq dx {x, Hp) = 1 because the rhs is fixed to 2J2N 
leptonic decay vr /U"^//^ for any factorization scale. 

The Sudakov functions can be written in terms of the momentum-dependent cusp anoma- 



lous dimension to read 17,18,19,20 



S (^,6,Q,/iF,/iR) 



1 /•MR=C2€Qd/i 



r 



cusp 



(7,c/(/^)) 



(7) 



where 7 = In (^^^7^) is the cusp angle, i.e., the emission angle of a soft gluon and the 
bent quark line after injecting at the cusp point the external (large) momentum by the 
off-mass-shell photon, and 



cusp 



(7,(?(/i)) = 21n 



A{g{fi)) + B{g{fi)) 



The functions A and B are defined by 



1 
2 



d 

cusp (gifi)) + /3ig)-^ICiCi, gifi)) 



2r 



TT 



+ -KC, 



TT 



and 



(8) 



(9) 



B (^7(/i)) = -- [IC {C,,g{fi)) + g (e, C,, gift))] 



7r 



(10) 



respectively, and the i^T-factor in the MS scheme is given by the expression 

10 



K 



18 6 



9 



■nfTp + /5oln(Cie^V2) 



(11) 



where Ca = = 3, Uf = 3, Tp = 1/2, 7e being the Euler-Mascheroni constant. The 
functions /C, Q are calculable within perturbative QCD |]T^|. Both anomalous dimensions, 
TcuspigifJ')) = Cf «s'^(/^^)/^ and 7q {g{^)) = —a^^{n'^)/7i, will be evaluated using the analytic 
model of |]I| in next-to-leading logarithmic order (for more technical details, we refer to [PT[|). 
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Fig. 1. Spacelike pion form factor calculated with cp^^ within different factorization schemes, except for 



the long-dashed line which shows the soft contribution |24|, in comparison with the existing experimental 
data |2^j2^. The lower solid line shows the IR-enhanced result obtained in our analysis, whereas the upper 
one stands for the sum of the long-dashed and the lower solid line. The dot-dashed line reproduces the 
calculation of |28|, which does not include an intrinsic /c_|_-dependence. 
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For simplicity, we follow |2T| and model the distribution of primordial (intrinsic) trans- 



verse momentum in the pion wave function in the form of a Gaussian normalized to unity: 

vl/^ (x, k^) = <p{x) P'g{x) exp [-g{x)(3'kl\ , (12) 

where g{x) = 1/xx and the quark masses are neglected. For 4>s.s{x) = 6xx, one has = 
0.883 [GeV~^] which corresponds to < >^/^= 350 MeV. 

Before we proceed with the presentation of our results, exposed in Fig. 0, let us at this 
point interject some comments regarding the role of the scales entering the calculation of 
the pion form factor. Whenever ^ < all Sudakov exponential factors are set equal 



to unity For all values of 6, there is a hierarchy of scales according to Aqqj-, ^ /if = 
Ci/b < fi^ = 6*2^(5 < Q- The limit /ir ~ fip can be interpreted as the minimum virtuality 
scale of exchanged quanta (or equivalently as the maximum transverse separation) in Th 
below which propagators cannot be treated within perturbation theory and are therefore 
absorbed into (f)-,^- In the present analysis we use /ip = Ci/b, /iR = C2^Q with ^ = x,x, y, y, 



and Ci = exp — |(27e — 1) 
(cf. [|I3|). For this choice o 



, C2 = l/\/2, the latter corresponding to the value C^"^ = 1 
the scheme constants, the logarithmic term in the K factor 
is eliminated. Note that for technical reasons, the argument of ag in Th is taken to be 

A^R = C2\f^Q, where ^ is either x or y. 

In order to show how the contribution to the pion form factor is accumulated in b space, 
we show in Fig. ^ the dependence of the scaled pion form factor against 6Aqcd- One 
observes a fast rise of the displayed curves as Q increases. Indeed, already for the smallest 
value shown, Q = 2 GeV, the form factor accumulates half of the whole contribution in the 
region b ~ 0.5/Aqcd- For still larger Q values, the form factor levels off already around 
b = 0.3/Aqcd for g = 5 GeV, and b = 0.25/Aqcd for Q = 10 GeV. This behavior of the 
curves uncovers how the IR stability and self-consistency of the perturbative treatment in the 



present scheme is improved, as compared to previous, conventional, approaches [Il4| , |21| , |28 . 

A similar expression to Eq. (Q) holds also for the main difference being that the 

latter contains only one pion wave function, and furthermore the associated short- distance 
part, Th, does not depend on ag in LO. The only dependence on the strong coupling constant 
at leading order enters through the anomalous dimensions of the cusp and the quark wave 
function. The result of this calculation is displayed in Fig. ^. Notice that in this case, we 
use /iR = C2^Q, where ^ = x or ^ = x. 

In summary, we have shown that modifying ctg in the IR region by a nonperturbative 
power correction, which removes the unphysical Landau singularity, may play a crucial role 
in the practical calculation of exclusive processes because it improves the IR stability of com- 
puted observables based on perturbation expansions without introducing external parame- 
ters to "freeze" the running strong coupling. Furthermore, we have given quantitative evi- 
dence that in this way it is possible to get an enhanced hard contribution to F-j^lQ"^), relying 
exclusively on the asymptotic form of the pion distribution amplitude, so that, though this 
contribution comprises Sudakov corrections and a primordial fc_L-dependence, it is not sup- 
pressed. Together with the soft part of F.„{Q'^), this contribution can account for the trend 
of the existing (admittedly low-accuracy) data without employing endpoint-concentrated 
pion distribution amplitudes. The same treatment yields for F^o^*^ a theoretical prediction 
which is in good agreement with the data. 
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Fig. 2. Dependence on the (critical) transverse distance b of the hard contribution to the scaled pion 
form factor, calculated with (^as within our IR-finite factorization scheme, for three different values of the 
momentum transfer: Qi = 2 GeV (solid hue), Q2 = ^ GeV (dashed line), and = 10 GeV (dashed-dotted 
line). 
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Fig. 3. Pion- photon transition form factor calculated with (/)as and IR enhancement (lower solid line). 



The other solid line shows the asymptotic behavior. The data are taken from |25 2' 
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